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Abstract 



We investigate static, spherically symmetric solutions of an Einstein- Yang-Mills-Chern-Simons system 
with negative cosmological constant, for an SO(6) gauge group. For a particular value of the Chern- 
Simons coefficient, this model can be viewed as a truncation of the five-dimensional maximal gauged 
supergravity and we expect that the basic properties of the solutions in the full model to persist in this 
truncation. Both globally regular, particle-like solutions and black holes are considered. In contrast 
Oh' with the Abelian case, the contribution of the Chern-Simons term is nontrivial already in the static, 

spherically symmetric limit. We find two types of solutions: the generic configurations whose magnetic 
gauge field does not vanish fast enough at infinity (although the spacetime is asymptotically AdS), whose 
mass function is divergent, and the special configurations, whose existence depends on the Chern-Simons 
£>C) . term, which are endowed with finite mass. In the case of the generic configurations, we argue that the 

divergent mass implies a nonvanishing trace for the stress tensor of the dual d — 4 theory. 

CO 

in 
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1 Introduction 

' It was originally found in d = 4 spacetime dimensions [I] , [2] , that a variety of well known features of asymp- 

totically flat self-gravitating non- Abelian solutions are not shared by their anti-de Sitter (AdS) counterparts. 
In the presence of a negative cosmological constant A < 0, the Einstein- Yang-Mills (EYM) theory possesses 
a continuous spectrum of regular and black hole non- Abelian solutions in terms of the adjustable parameters 
that specify the initial conditions at the origin or at the event horizon, rather than at discrete values of 
these parameters. The gauge field of generic solutions does not vanish asymptotically, resulting in a nonzero 
magnetic flux at infinity. Moreover, in contrast with the A = case, some of the AdS configurations are 
stable against linear perturbations [3]. 

As found in [3], [5], some of these features are shared by higher dimensional EYM solutions with AdS 
asymptotics. Since gauged supergravity theories generically contain non-Abelian matter fields in the bulk, 
these configurations are relevant in an AdS/CFT context, offering the possibility of studying some aspects 
of the nonperturbative structure of a CFT in a background gauge field (BJ. On the CFT side, the boundary 
non-Abelian fields correspond to external source currents coupled to various operators. 

Given its relevance in the conjectured AdS/CFT correspondence [7], jS], the case of Af — 8, d = 5 gauged 
supergravity [S], [TUj is of particular interest. The bosonic sector of this theory consists of the metric, twenty 
scalars and fifteen 50(6) Yang-Mills gauge fieldfl Apart from the usual F 2 term, the Yang- Mills (YM) 
fields have in this case a non-Abelian Chern-Simons (CS) term in the action, which unlike in the Abelian 
case does not vanish when subjected to spherical symmetry. 



1 Note that the field content of the full Af = 8, d = 5 gauged supergravity is richer. However, a number of bosonic fields can 
be consistently set to zero [101 . 
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Solutions of this model have been considered by several authors for various consistent truncations, with 
subgroups of 50(6) (see e.g. the recent work [TT] and the references therein). However, to our knowledge, 
no attempt has been made to construct non-Abelian solutions for the general case of the full 5*0(6) gauge 
group. In particular, the effects resulting from the introduction of the CS term have so far not been studied. 

This paper is aimed as a first step in this direction, by taking a truncation of the Af — S, d = 5 
model corresponding to a pure Einstein- Yang-Mills-Chern-Simons (EYMCS) theory (i.e. with a negative 
cosmological constant, but with no scalar fields). We propose an Ansatz for a spherically symmetric 50(6) 
gauge group and investigate the basic properties of both the black hole, and, particle-like globally regular 
solutions. Special attention is paid to the new features induced by the CS term. 

As originally found in [12], [4], [5], a generic property of higher dimensional EYM solutions is that 
their masses and actions, as defined in the usual way, diverge. (For a recent review of these solutions, see 
[13].) This can be understood heuristically by noting that the Derrick scaling requirement is not fulfilled in 
spacetimes for dimension five and higher. To our knowledge, the only mechanism for regularising the mass of 
the d > 4 asymptotically flat or (A)dS non-Abelian gravitating solutions, proposed so far in the literature, is 
to include higher order terms in the YM hierarchy [M], [5], [15] and the corresponding YM-Higgs terms 16JH1- 
These are the YM counterparts of the Lovelock gravities (or the hierarchy of Einstein systems) , and occur 
in the low energy effective action of string theory [TBI [19] . 

One of the main features of the present work is the introduction of the CS term in 4 + 1 dimensions, as an 
alternative to the higher order curvature terms of the YM hierarchy employed previously for regularising the 
mass. It turns out that this prescription does result in finite mass solutions, but in addition to these, we find 
solutions with divergent mass in what we have termed as the generic case. These last are characterised by 
a continuum of values of the shooting parameters, which is a typical feature of EYM system with a negative 
cosmological constant. The finite mass solutions on the other hand, termed as the special case, are special 
in that they occur only for a discrete set of values of the shooting parameters. Of course, in the absence of 
the CS term the only solutions that exist are ones with divergent mass. 

Although the spacetime still approaches asymptotically the maximally symmetric AdS background, the 
mass and the total action of a generic solution present a logarithmically divergent part. The coefficient of 
the divergent term is fixed by the square of the induced non-Abelian fields on the boundary at infinity. We 
shall argue that the logarithmic divergence of the non-Abelian AdSs configurations does not signal a problem 
with these solutions, but rather provides a consistency check of the AdS/CFT conjecture, the coefficient of 
the divergent term in the action being related in this case to the trace anomaly of the dual CFT defined in 
a background non-Abelian field. Moreover, one can define a mass and action for the generic solutions by 
using the counterterm prescription of [20] . The counterterms here depend not only on the boundary metric 
but also on the induced non-Abelian fields on the boundary. 

However, perhaps the most interesing feature of the EYMCS model is the existence of a set of solutions 
with finite mass. In the case of these solutions, as for the well known d — 4 Bartnick-McKinnon solitons [21., 
they exist only for discrete values of the shooting parameters (associated with the initial values of the gauge 
fields). As can be seen by using a simple Derrick- type argument, they are supported by the contribution of 
the non-Abelian CS term, a prescription which can be exploited only in odd dimensional spacetimes where 
a CS term is defined. 

The paper is structured as follows: in Section 2 we present the general framework and analyse the field 
equations and boundary conditions. We present the numerical results in Section 3, special atention being 
paid to solutions with a finite mass. The computation of the mass and electric charge of the solutions is 
addressed in Section 4. We conclude with Section 5 where the significance of, and further consequences 
arising from, the solutions we have constructed are briefly discussed. 

2 It is in principle possible to supply such higher scaling terms by employing only Higgs kinetic terms, or, the kinetic terms 
of suitably gauged higher dimensional sigma models I17| . but these have not been attempted. 
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2 The model 



2.1 The action 

We consider the following action 



d*xV=9 ( tt^(R - 2A) - Cym ) - I d b x C cs , (1) 

M 



M \16irG 



where 



Cy M = ilkiV"}, (2) 

is the usual Yang-Mills lagrangian for a gauge group 50(6) (with = d^A u — d u A^ + eL4 M , A v \ the gauge 
field strength tensor), and 



pnuppa _ eF ^A p A' 7 + -e 2 A p A u A p A a 
5 



(3) 



is the CS ternjf] (with n the CS coupling constant), while A = —6/£ 2 is the cosmological constant and e is 
the gauge coupling constant. 

These are basic the pieces which enters the bosonic action of the d = 5, Af = 8 gauged supergravity [S], 
[TU] . the CS coefficient being k = 1/8 in this case. However, the full N — 8 system contains in addition 
twenty scalars, which are represented by a symmetric unimodular tensor. These scalars have a nontrivial 
potential approaching a constant negative value at infinity which fixes the value of the effective cosmological 
constant. Although ignoring the scalar sector is not a consistent trucation of the general Af — 8 model, we 
expect that the basic properties of our solutions hold also in that cas^|. 

The field equations are obtained by varying the action ([T]) with respect to the field variables , A^ 

R^-^g^R + ^g^ = SnGT^, (4) 



v-g 

where the energy momentum tensor is defined by 

Tfi V - Tt^F^FvpgP* - X - 9ilv F a pF a ?y (5) 
One can show that this tensor is covariantly conserved (i.e. V AI T'" y = 0) for solutions of the YMCS equations. 
2.2 The spherically symmetric Ansatz 

In this work we shall restrict to simplest case of static, spherically symmetric solutions. Thus we consider a 
metric Ansatz in terms of two metric functions N(r) and cr(r) 

ds2 = W{r) + r2dn " ~ N{r)a2 {r)dt2 ' (6) 



where we have found convenient to define 



AW-i-=gi + £. (7) 



3 The factor of i appears in Q because we are using an antihermitian representation for the SO (6) algebra matrices. 

4 This is the situation for d = 4 EYM solutions. As discussed e.g. in [22] . the properties of the EYM-dilaton solutions (the 
dilaton field possessing a nontrivial potential approaching a constant negative value at infinity) are quite similar to those of the 
pure EYM-AdS case. 
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the function m(r) being related to the local mass-energy density (as defined in the standard way) up to some 
factor, r and t are the radial and time coordinates, while <if2§ is the metric on the round three-sphere. 

The static, spherically symmetric SO(6) YM fields are taken in one of the two chiral representations of 
SO(6), such that the spherically symmetric Ansatz is expressed in terms of the representation matrices, 



(8) 



where Sj = — Sj = iji , £5 = — £5 = 275 , T.q = +^6 = 2, are defined in terms of the usual Dirac gamma 
matrices 7.;. 

Our spherically symmetric Ansatz for the SO(6) YM connection = (A t ,Ai) is a variant of Witten's 
Ansatz for the axially symmetric instanton [24]. Also, this Ansatz is one of the two SU(4) Ansatze proposed 
in [25], namely the one employing Dirac gamma matrices as opposed the one employing the Gell-Mann 
matrices^. It is expressed as 



At = - e {-(e X (r)) M x 3 ^ M -x 7 (r)^ 56 ) 



(9) 



f 



7 (r) + 1 



T, jM + A 7 r (r) Xi S 5 6 



In the above relations 



where Xi = Xi/r (with x 1 the usual Cartesian coordinates on R and XiXi = 1 
i,j = 1, 2, 3, 4 and the index M runs over 5, 6. Also, e is the two dimensional Levi-Civita symbol. 

After taking the traces over the spin matrices, it is convenient to relabel the triplets of radial function 



as 



),X= (x M ,X 3 ) and A r = (A^,A 3 r ), with M = 1,2 now. 



Substituting © , (flTJj) in the YM Lagrangian density we find the compact expression for the reduced one 
dimensional reduced YM action densitMj 



-ra I V D, 



1 



1 r J 



(10) 



The calculation of the reduced CS action density is rather more tedious since unlike (fTU|) . this term is not 
gauge invariant. It can be expressed in a compact way as 



-cs 



4c. :! 



12(\4> d \ 2 - l)A a {e abc x b 4> c ) + 3|/| 2 {4> a DrX a ) 



+ 6(0V)(0 a ^ a )-9|0 b | 2 (x a ^0 Q ) 
- (24> 3 + 5WD rX a ) + (3|^ a | 2 - 2i? - 5) D rX 3 

+ (6x'V + 7 X a ~ 20V) D 7 4 a - 2(x a r + x 3 ) D r <p 3 



(11) 



which of course does not feature the metric functions. Note that (fTTj) is not a scalar after contraction of the 
indices (a, 6, c). This is a consequence of the gauge variance of the CS density. 
In both ([TO]) and (fTTj) we have used the notation 



D r (f> a = d r cj) a + e abc A'. 



D rX a = d rX a +e abc A h r x c , 



(12) 



which are SO(3) covariant derivatives of the two triplets 4> = cf) a — (</> M ,</> 3 ), and x = X a — (x M iX 3 )i with 
respect to the S*0(3) gauge connection A r = A". But A r is really a pure — gauge since in one dimension there 
is no curvature. As such, it can be consistently set equal to zero. But more importantly, taking the variations 
5A r leads to the constraint equations, which are first integrals of the equations for <p and x, an d which play 



5 This distinction is important since the Dirac gamma matrix Ansatz cannot be contracted to a ST/ (2) subalgebra, while 
clearly Gell-Mann matrix Ansatz does have a 5(7(2) subalgebra. 

6 Here we ignore the (irrelevant) angular part in y/—g. 
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an important technical role in the numerical integrations. We will return to these below. The ocurrence of 
constraint equations in a system supporting what are basically sphaleron solutions is completely expected, 
as the solutions we construct are indeed sphalerons, just like the familiar Bartnik-McKinnon solutions. 
Needless to say, the consistency of the Ansatz used has been verified, so it is sufficient to work with the 
reduced one dimensional Lagrangian (fTU]) - (fTTj) . 

Finding solutions within the general YM Ansatz ©, which after setting A r — still features six in- 
dependent functions, is technically a difficult task. A further consistent trucation of the general Ansatz is 
(f) 2 — x 2 — 0, leading to an EYMCS system with six unknown functions, four of them being gauge potentials 
parametrising the gauge field, and, two metric functions. Indeed, the two gauge functions suppressed are 
redundent and would only be excited in an eventual stability analysis of our sphalerons. 

To make connection with notations used in previous work [H [12j [14] on d = 5 EYM solutions, we adopt 
the notation 

<f> 1 (r)=w(r), ci> 3 (r)=w(r), x 1 (r) = V(r), x 3 (r) = V(r). (13) 

The resulting system has some residual symmetry under a rotation of the 'doublets' w(r),w(r) and 
V(r),V(r) with the same constant angle u {e.g. w — > wcosu + wshiu etc.) One can use this symmetry 
to consistently set w(r) = V(r) = (or w(r) = V(r) = 0) which results in a particular truncation of 
the system, which we shall exploit in Section 3.2. Note that for configurations with w(r) — V(r) = 
the gauge potentials are invariant under the "crural" transformations generated by £5. The configurations 
with w(r) = V(r) = instead change just by a sign under the same transformations. Also, this Ansatz is 
invariant under the parity reflections transformation <fi a — > —<fi a , x a ~~ * ~X a ■ The asymmetry beween (w, V) 
and (w, V) is manifested by the different set of boundary conditions they satisfy. 



2.3 The equations and boundary conditions 

Inserting this Ansatz into the action ([T]), the EYMCS field equations (j4|) reduce to (to simplify the notation 
we denote a 2 = 167rG/(3e 2 ) and absorb a factor of 1/e in the expression of n): 

m> = \a 2 (3r( N ( W > 2 + uP) + ^ + g ' ^ ) + J (V 2 + V' 2 + ^(Vw - Vwf 

a' 3a 2 f ,o _,o 1 

— = w' 2 + w' 2 H - 

a 2r \ N 2 a 

'2w(w 2 + w 2 - 1) V 



w ' 2 + ™' 2 + ~ V ™) 2 ) ' ( 14 ) 



(raNw')' = ra \^ W[W ±^ '- + -^(Vw - Vw)J + 4k (v'(w 2 + w 2 - 1) + 2w'(Vw - Vw)^j 

(2w(w 2 + w 2 — 1) V \ f~ ~ \ 
r 2 + 'a 2 N <yVW ~ V ^ ) + 4K V ^ + ^ ~ ^ + 2w '( Vw ~ V{b " i ) 

w(Vw - Vw) + 12k(w 2 + w 2 - l)w\ 
w(Vw - Vw) + 12k(w 2 + w 2 - l)w', 



- a ' aN 
,r 3 V\, 3r 



(7 aN 
together with the constraint equation 



r 



3 



(VV - VV) + 3rNa(ww' - ww') - 12k(Vw - Vw)(w 2 + w 2 - 1) = , (15) 



<7 

which originates from the variational equation for 8A r (where a prime denotes a derivative with respect to 
r). 

These equations support both globally regular and black hole solutions. The only known closed form 
solutions of these equations are discussed in the next subsection and are trivial in some sense, since the 
magnetic gauge potentials do not feature any dependence on r. However, it is concievable that non-Abelian 
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analytic solutions can be found by studying the first order Bogomol'nyi equations of the full Af = 8 gauged 
supergravity model, with all scalar functions included. This was the case of other gauged supergravity 
theories, the most famous example being the Chamseddine-Volkov solution [33] of the Af = 4, d = 4 
Freedman-Schwarz model [27]. One might therefore expect the full Af = 8 model to support BPS solutions 
describing also non-Abelian globally regular solitons, which actually constrasts with the case of an Abelian 
truncation. However, given the large number of matter functions, even finding the explicit form of the first 
order Bogomol'nyi equations of the Af — 8 supergravity model with non-Abelian fields is bound to be a 
very difficult task, which has not been addressed so far in the literature. This is not surprising since the 
analogous task in the construction of the Chamseddine-Volkov solution involves simply a 577(2) gauge field 
and a single dilaton field. 

However, one can analyse the properties of the solutions of the system (114|) by using a combination of 
analytic and numerical methods, which is sufficient for most purposes. 

The globally regular configurations are nontrivial deformations of the AdSs and have the following ex- 
pansions near the origin r = 0: 

/■ \ -, .2 ^/ 4x ~ / \ gi(V(0) — 24nba ) , - 
w(r) = 1 - br 2 + 0(r 4 ), w(r) = * y v ; 2 —r A + 0(r 5 ), 

V(r) = V(0) +6b 2 n<7 r 2 + 0(r 4 ), V(r) = gi r + 0(r 3 ), (16) 



a 



2 ( 5l 2 +66V 2 ) , , 6 , , 3a 2 (g 2 + 4b 2 a 2 ) , 



= " vai -u/ r 4 + G(rb)j cr(r) = (To + w + ^r 2 + 0(r 4 ). 

2cr Q 4cr 

The free parameters are b — —^w"(0), V(0), gi = V'(0) and cto = cr(0). The coefficients of all higher order 
terms in the r — ¥ expansion are fixed by these parameters. 

We are also interested in solutions having a regular event horizon at r = > and representing 
non-Abelian generalisations of the Reissner-Nordstrbm-AdSs (RNAdS) black hole. To simplify the general 
picture we shall consider mainly nonextremal black holes, in which case N(r) has a single zero at r = and 
&{ r h) > 0. We expect that extremal black holes also exist for the full 50(6) theory, but we have restricted 
their numerical construction only to the particular truncation (w = V = 0) of the system, alluded to in 
the previous subsection, which is of course, a consistent truncation. We have made this restriction simply 
due to our desire to render the numerical task easier. From our study of this particular truncation of the 
system, we deduce that it is likely extremal black holes exist also for the full 5(3(6) theory. 

For the nonextremal case, the field equations imply the following behaviour as r — > r^ in terms of five 
parameters w h = w(r h ),w h = w{rh), V h = V(r h ), Vi = V'{r h ) and a h = a{r h ): 

w(r) = w h +Wi(r — r h ) + 0(r - r h ) 2 , w(r) = w h + w\{r - r h ) + 0(r - r h ) 2 , 

V(r) = V h + V 1 (r-r h ) + 0(r-r h ) 2 , V(r) = ^ + ^(r - r h ) + 0(r ~ r h ) 2 , (17) 

2 n , r 2 h a 2 f rlV 2 (w 2 + w 2 ) 3(w 2 + w 2 - if 

cr(r) = o h +oi(r - r h ) + 0(r - r h ) 2 , 



m{r) = rf L (l + + - h 1 \> l ' > 1 ' + y > l ' > l '- ) (r - r h ) + 0{r - r h )\ 



where 



Wl 



Ar h £ 2 a h wl(2Kr h Vi + a h w h ){w 2 h + w 2 h - 1) 



<7l = 



-Ar 2 (2r 2 + £ 2 )a 2 w 2 + a 2 1 2 (3o 2 w 2 ( W 2 + w 2 I) 2 + r{V 2 {w 2 + w 2 )) ' 

- 4r h £ 2 OhW h w h (2Kr h Vi + (JhW h ){w 2 h +w 2 h -l) 

Wl ~ -K^K + ^ 2 e 2 <J 2 w 2 (w 2 + w 2 l)^ + rK-Salw 2 + a 2 £ 2 V 2 (w 2 + w 2 )) ' 
2^a 2 r h l i (rlw 2 {2nrhVi + a h w h ) 2 (w 2 + w 2 - l) 2 (w 2 + w 2 ) 



> h uj h - 3a 2 £ 2 a 



Whi< +<- I) 2 + - a 2 PV 2 {w 2 h + w 2 ))^ ' 
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2.3.1 Large r asymptotic expansions 

The expansion at infinity of the solutions is more complicated and involves separate analyses in the generic 
and the special cases. 

In the generic case, the potentials parametrising the non-Abelian gauge field take arbitrary values as 
r — > oo, with the leading order behaviour 



w 2 . logr , _ , U)2 , _ log 



r 



w(r) = w + — + w c — — + . . . , w(r) = w + — + w c —^- + . . . , 
V ( r ) = y o + ^ + ^^¥ : + ---' V{t) = Vo + ± + V^ + ..., (19) 



with 



i(r) = Mo + \a 2 ((tag + fig - l) 2 + £ 2 (V w - V w ) 2 ) logr, a(r) = 1 - 

£ 2 i ~ \ 

wc = ~ [2w (w 2 +w% - 1) + {V Wo - VqWq)J 



a 2 (w 2 + w 2 c ) log 2 r 



w c = -- (2w (wg + tSg - 1) + Fo^ 2 (Vb^o - Vq* q )J , (20) 

3 3 
V c = --t w (V w - V w ), V c = --e 2 w (V wo - V w ), 

where Wq,Wq,Vq,Vq and W2,W2,q,q are arbitrary parameters satisfying the constraint 

3(w 2 w - w 2 w ) + e 2 (qV - qV ) + 6£ 2 k(V q w - V w )(w 2 + tog - 1) = 0. (21) 

Thus, similar to the well known case of a SO(3) gauge group 0], the generic non-Abelian configurations 
have a nonvanishing magnetic field on the AdS boundary (i.e. F^F^^^^ ^ 0). As a result, one can see 
from the above relations that the mass function m(r), and hence also the action of these solutions, diverge 
logarithmically Q- 

However, despite the divergence of the mass, the spacetime is still asymptotically AdS, the large r 
behaviour of the metric function N(r) being N(r) — > r 2 /£ 2 + l. Asymptotically AdS solutions with diverging 
mass have been considered recently by some authors, mainly for a scalar field in the bulk (see e.g. [311). 
In this case it might be possible to relax the standard asymptotic conditions without loosing the original 
symmetries, but modifying the charges in order to take into account the presence of matter fields. In Section 
5 of this work we shall argue that this is also the case of the EYMCS configurations with the general 
asymptotics (|19p . By using a counterterm approach, one can define a mass for these solutions, which is 
fixed by the parameter Mq appearing in (119[) . (Note that for generic solutions not only ro(r) diverges 
logarithmically as r — > oo but also the terms r 3 V'(r) and r 3 V'(r) which, as argued in Section 4, fixes the 
electric charge(s) of the solutions. In the numerics, we have studied mainly the solutions with V c = V c = 0.) 

For the special configurations, which support finite mass, the required asymptotic behaviour at large r 
is \<Pq \ — > 1 and \<j> x x\ — > (i.e. w" 2 +Wq — > 1 and WqVo — WoVo — > 0). These conditions can be satisfied only 
in the presence of the CS term. Different from the case of a simple EYM theory [3] , the existence of finite 
mass configurations here is not forbidden by the Derrick-type scaling argument. Indeed, the numerics in the 
following Section indicate the existence of a subset solutions with finite mass, with the following expansion 



7 The existence of a logarithmic divergence in the action is a known property of some classes of AdSs solutions which are 
endowed with with a special boundary geometry 28_ . The coefficients of the divergent terms there are related to the conformal 
Weyl anomaly in the dual theory 29 30 . However, this is not the case for the non-Abelian AdSs configurations here, which 
have the same boundary metric as the Schwarzschild-AdS solution and hence feature no Weyl anomaly in the dual CFT. 
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at infinity 

w(r) = since + ^ + 0(l/r 4 ), u)(r) = cosci + ^ + 0(l/r 4 ), 

l/(r) = $sina + -| + 0(l/r 4 ), F(r) = $cosa + 4 + 0(l/r 4 ), (22) 
m(r) = M - *V + ^ 2 +^f + ^)) + ^ ^ = x _ a^f^) + (l/r«), 

where the amplitude of the electric potential at infinity is fixed by 

3(w2 cos a — W2 sin a) . . 

* = -^7 — f. (23) 

t^yq cos a — q sin a) 

Thus the free parameters in the far field expansion of this special set of solutions are Mq, a = 
arctan(y(co)/V r ((xi)) and the coefficients q, q, W2,W2 of the 1/r 2 decaying terms in the non-Abelian po- 
tentials. 



2.4 Particular cases 

The simplest solution of the field equations has pure gauge fields {F^ u = 0) and corresponds to the 
Schwarzschild-AdSs black hole 

r 2 M 

iV(r) = 1 + J ' cr ( r ) = 1 ' w(r) — sin a, w(r)=cosa, F(r)=$sina, F(r)=$cosa, (24) 

where $, a are arbitrary constant. 

The embedding of the RNAdS Abelian solution is recovered for 

v 2 Ad oi 2 a 2 ~ a 

AT( r ) = l+--— + -f- , a(r)=l, w(r)=V(r)=0, w(r) = ±1, V(r) = $ + ±. (25) 

The only AdS exact solutions with nontrivial non-Abelian fields known so far is: 
r 2 Mn + — log r 

N{r) =1 + ^2 " — -i — > CT ( r ) = 1 ' W W = *W = °> y W = y W = °> ( 26 ) 

(with Mo an arbitrary positive constant). This solution was obtained in [3] and describes a Reissner- 
Nordstrom type geometry in EYM theory with a gauge group 50(3) (note that the mass function is loga- 
rithmically divergent in this case). Its embedding in the 50(6) gauged supergravity model and the extremal 
limit has been discussed in the recent work 

A particularly interesting model is found by taking w(r) = V(r) = (or equivalently w(r) = V(r) = 
0). This is our particular truncation of the full 50(6) model. The resulting solutions are those of the 
SU(2) x U (1) truncation of the modeH parametrised in terms of the representations of the algebra of 5[/(4) 
instead of 50(6). (Of course in that case both gauge groups have the same gauge coupling constant.) One 
can see that the CS term is still nontrivial in this case 

C cs = i KV(r)e aia2a3ai Tr {F a ^F a ^} = 12 nV (r)w' (r)(w 2 (r) - 1), (27) 

(with dj = 1, . . .4). As we shall argue in the next Section, the solutions of this particular truncation contain 
already the basic features of the full model. However, they are much easier to study numerically. 
The YM-CS equations in this case admit the first integral 

V'(r) = ^(K + 4kw{w 2 -3)), (28) 
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' v ' V(r)/5 

-2-1012 _1 ' 5 -0.5 0.5 1 1.5 2 2.5 

Figure 1: The profiles of a typical globally regular solution (left) and a black hole solution (right) of the EYMCS 
equations are presented as a function of the radial coordinate r. For these solutions, the mass function m(r) diverges 
logarithmically as r — > do. 



with K being an integration constant. One can easily see that the solutions are regular at the origin, r = 0, 
only if K = 8k. The value of K is not fixed a priori for black hole solutions. 

The asymptotics of the SU(2) x U(l) solutions can easily be read from the general relations (fTo]) , (IT91 
and (|22l) . At infinity, the finite mass solutions have a = ±ir/2 in (|22|) . with the gauge potentials 

F ( r ) = *-^^ + °( 1 A 4 ) I w(r) = ±l + ^ + 0(l/r 4 ). (29) 

Thus, for the case w(oo) = —1 studied in this work, the parameter q fixing the Abelian electric charge of 
these solutions is q — — (4k + K/2) for black holes and q = —8k for globally regular solutions. 



3 Numerical results 

We start by noticing that the equations (Q~4|) are not affected by the transformation: 

r -> Xr, m -> A 2 m, I -> XI, V -> V/A, V" -> 7/A, a -> Aa , k -> Ak, (30) 

while w, w and a remain unchanged. It follows that one can always take an arbitrary positive value for a. 
The usual choice is a = 1, which fixes the EYM length scale L = -\/87rG/(3e 2 ), while the mass scale is fixed 
by Ai — 87r/(3e 2 ). All other quantities get multiplied with suitable factors of L. However, in this Section, 
to avoid cluttering our expressions with a complicated dependence on (G, e), we fix the value of a at a = 1, 
and ignore the extra factors of e and G in the expressions of various global quantities. 

Therefore the remaining input parameters are the AdS length scale t and the CS coupling constant k. 
Determining the pattern of the solutions in the parameter space represents a very complex task which is 
outside the scope of this paper. Instead, we analysed in detail a few particular classes of solutions which, 
hopefully, reflect all relevant properties of the general pattern. For dcfinitcness we set I = 1 in our numerical 
analysis, although we have found nontrivial solutions also for other values of the cosmological constant^. 

8 This truncation of the full SO (6) model shares a number of common features with the five dimensional N = 4 gauged 
SU(2) X U(l) supergravity model considered in 1321 . For example, a first integral similar to H28II appears there also. However, 
the solutions in 1321 have an extra dilaton field with a Liouville potential and thus are not asymptotically AdS. 

9 In particular, the finite energy solutions survive in the limit A — > 0, being supported by the CS term (this contrasts strongly 
with the case of a pure EYM theory 1121 ). A discussion of the asymptotically flat EYMCS solutions will be presented elsewhere. 
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Figure 2: A number of relevant parameters are plotted as functions of the coefficient b for globally regular solutions 
of the EYMCS model with £ = 1, k = 1. These solutions have finite electric charges, which are fixed by q and q. 



Instead we have looked for the dependence of the solutions on the value of the CS coefficient k, which has 
not been fixed a priori. (This investigation has been partly motivated by the study in |33j of the Einstein- 
Maxwell-CS system, which revealed a nontrivial dependence of the properties of the solutions on the value 
of k. We shall see that this is also the case for the solutions constructed in this work, which feature a critical 
value of CS coefficient.) 

The resulting set of six ordinary differential equations^ is solved with suitable boundary conditions 
which result from (fTrJ|) . (fTT|) , (|19p and (j2"2"|) . The numerics employs a collocation method for boundary- 
value ordinary differential equations equipped with an adaptive mesh selection procedure |34j . Typical mesh 
sizes include 10 3 — 10 4 points. The solutions have a relative accuracy of 10~ 7 . In addition to employing 
this algorithm, some solutions were also constructed by using a standard Runge-Kutta ordinary differential 
equation solver. In this approach we evaluate the initial conditions at r = 10 -5 (or r = + 10~ 5 ), for global 
tolerance 10~ 12 , adjusting for shooting parameters and integrating towards r — > oo. We have confirmed that 
there is good agreement between the results obtained with these two different methods. 

The properties of the solutions depend on the input parameters, but it is rather difficult to find a general 
pattern. However, a feature shared by all asymptotically AdS solutions is that the metric functions m(r), 
cr(r) monotonically approach their asymptotic values, which can easily be seen from the corresponding 
field equations. Also, so far we could not find solutions where the electric potentials V(r), V(r) present 
oscillations, even though such solutions are allowed. 



3.1 The generic solutions 

The generic solutions studied here are for the d = 4 + 1 dimensional model, which features a Chern-Simons 
(CS) term. But solutions with similar properties are found also in the EYM model with no CS term. All 
these solutions bear a qualitative similarity to those of the more familiar d = 3 + 1 EYM model [1] , [2] . These 
solutions can also be seen as higher gauge group generalisations of the EYM solutions in [4], but unlike the 
latter they feature a nontrivial electric potential, made possible by the larger gauge group. (The electric 
potential necessarily vanishes for a d = 5 static, spherically symmetric SU(2) gauge field). 

Considering first the case of globally regular configurations, one finds that solutions approaching asymp- 
totically the AdSs background exist for compact intervals of the initial parameters w"(0), V(0), V'(0) and 
<t(0). The values of the parameters wo,wo,Vq,Vo, u>2, u>2 and q which enter the asymptotics of the solutions 

10 Although we have solved the second order YM equations in ||14|I . we have also monitored the constraint (H) , which was 
always satisfied with very good accuracy. Also, the equation Ill5jl has been used to construct the asymptotic expansions JIB), 
CE3, CE9} and (22). 



10 



branch 2 



0.1 0.15 0.2 0.25 0.3 0.1 0.15 0.2 0.25 

K K 

Figure 3: A number of relevant parameters are plotted as functions of the Chern-Simons coupling k for finite mass, 
globally regular solutions of the EYMCS model. One can see that new branches of solutions emerge as k increases. 



are fixed by the numerics. (There are also branches of solutions with a different asymptotic behaviour, 
which stop to exist for finite values of r. To study these configurations, one needs to employ a metric Ansatz 
different from ©. Such solutions, being not asymptotically AdS, are of no interest here.) 

Of the full set of solutions with AdS asymptotics, we have paid special attention to the physically more 
interesting case of configurations with a 1/r 2 decay of the electric potentials at infinity (i.e. V c = V c = and 
VqWq — VqWq = 0), this being the only case reported in this Section. These solutions have a finite electric 
charge, although their mass functions will diverge asymptotically since \(f>\ — \/u>q + Wq -/> 1 here. 

A typical configuration with a regular origin is presented in Figure 1 (left), for k = 1. One can see that 
the mass function diverges logarithmically while <r(r), w(r),w(r) and V(r),V(r) asymptotically approach 
some finite values. Solutions with nodes in w(r), w(r) were also found. 

In Figure 2 we plot a number of relevant parameters as a function of the coefficient b in the initial data 
at r — (b — — w"(0)/2), for a family of asymptotically AdS solutions. (One of the parameters there is 
Mq appearing in (j!9[) . which in Section 5 we argue that it can be taken as the renormalised mass of the 
solutions; note that Mq may take also negative values). This branch ends for some finite values of b, where 
~9tt(0) = c(0) — >• while Vo,q diverge. The condition \<f> x x\ — > as r — > oo has been enforced by treating 
V(0) as a shooting parameter. Then ^'(0) is a free parameter while <r(0) results from the numerics (the 
solutions in Figure 2 have V'(0)/(Jq = 0.15). 

The results in Figure 2 show that the for generic solutions w 2 (oo) + w 2 (oo) -ft 1. From (flU)) . this leads 
to a divergent mass-energy as defined in the usual way. However, one can see that the condition \cj)\ — > 1 is 
satisfied for a discrete set of the parameter b [e.g. b ~ 0.3219 x 10~ 3 and b ~ 0.6125 x 10~ 5 for the data in 
Figure 2). This suggests the existence of several branches of finite mass solutions parametrised by V^'(0) (or, 
equivalently, V(0)), which is confirmed by the results in the next subsection. 

Black hole solutions have been found as well, presenting the same general features. Here also we find a 
continuum of solutions with arbitrary values of gauge potentials at infinity, the relevant parameters being 
the values of the gauge potentials at the event horizon as given by (jTT)) . Again, finite mass black holes are 
found only for specific values of the gauge potentials on the horizon. 

As a general remark, we note that the presence of the CS term is not crucial for the existence of the 
generic solutions (i.e. with a divergent mass). We have found solutions with rather similar properties also 
for k = 0. Thus, the role of the CS term is indispensable only for the construction of special, finite mass, 
solutions to be presented in the next subsection. 
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Figure 4: The profiles of two typical globally regular EYMS solutions are presented as function of the radial 
coordinate r. m(r) and a(r) are metric functions, which w(r),w(r) and V(r),V(r) are non-Abelian potentials. 



3.2 Finite mass solutions 
3.2.1 Regular configurations 

In the numerics, special attention has been paid to solutions with a finite mass. As noted above, in this 
case, three of the four parameters in the data at r — are fixed. The remaining free parameter was chosen 
to be V(0) (similar results are found when V(0) is chosen instead). 

Remembering the invariance under the parity reflection cj> — > —0, x ~^ ~X it turns out to be sufficient to 
consider 7(0) > (or V(0) > 0). 

Thus, corresponding to a choice of the coupling constant k and of the cosmological constant A, there 
exist in principle a family of finite mass charged EYMCS solutions labeled by the value at the origin of one 
of the electric potentials, in this case, V(0). The pattern of these solutions turns out to be extremely rich, 
with some unexpected features. 

In order to illustrate this, we first fix V(0) = 0.3 and study the solutions as functions of the CS parameter 
k. (Qualitatively, the same results have been found when considering other values of V(0)). The numerical 
results show that a branch of solutions with w(r) = V(r) = always exists for sufficiently large values of k, 
for instance, k > ko ~ 0.07 in the present case. These are the solutions of the reduced St/ (2) x t/(l) model 
with the asymptotic angle a = ir/2. For convenience, we will refer to this branch as the main branch. In 
the limit k — > kq the metric function cx(r) vanishes at the origin and the solution becomes singular. For all 
values of the other parameters, no finite mass solutions have been found for k < kq. 

The interesting feature is that new branches of solutions with notrivial functions w(r), V(r) emerge from 
the main branch at critical values of k. In our cases, the first branch of excited solution appears n ~ 0.0953 
and a second branch at k ~ 0.1875. This behaviour is illustrated in Figure 3, where a number of relavant 
global parameters are plotted a functions of the CS coupling constant k. It should be noted that, for a fixed 
k, the excited solutions have larger masses than the solutions on the main branch. 

The profiles of the two first excited solutions corresponding to k = 0.2 is presented in Figure 4. One can 
sec that the functions w, w parametrising the magnetic field of the excited solutions develop more pronounced 
oscillations before becoming constant in the asymptotic region. 

It is also natural to study the spectrum of solutions in terms of one of the charges, say q for a fixed 
value of k. The second electric charge q is determined from numerics. The results of our analysis for the 
value k = 0.2 are summarised in Figure 5. In this case, two excited solutions are available. Several relevant 
parameters are plotted there versus q (see Eqn. <\22^i ) for two excited solutions. Fixing the parity symmetry 
by means of ^'(0) > 0, the numerical analysis reveals that the solutions of the branch "1" (respectively 
"2") are characterised by negative (respectively positive) values of q. Also, they exist up to a minimal 
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Figure 5: A number of relevant parameters are plotted as functions of the electric charge q for finite mass, globally 
regular solutions of the EYMCS model with a Chern-Simons coefficient k — 0.2. The dashed and solid curves denote 
different branches of solutions. 



(respectively maximal) value, say q = q cr . For k = 0.2 we have q cr s=s —7.8 and q cr « 2.5 respectively, for 
the first and second branches. In the limit q cr — > 0, the excited solutions converge to the main solution. The 
ending of the branches at q = q cr is more subtle. Indeed, our numerical results show that another family 
of solutions (with larger mass) exists in the region \q\ < \q cr \, backbending from the branch coming directly 
from the main solution. These new branches are shown on Figure 5; however, we have not attempted to 
construct further branches in this region, although they are likely exist. 



3.2.2 Black holes 

The EYMCS system presents also black hole solutions which were constructed using similar techniques. In 
contrast with the regular solutions presented above, the main thrust here is confined to the solutions of 
our particular truncation of the full SO (6) model. This restriction is made to simplify an otherwise very 
complex numerical task. 

The Hawking temperature and the entropy of the black holes are given by 

a(r h )N'(r h ) A H 3 

T H = — , S = —, with A H = V 3 rl, (31) 

where V3 = 2ir 2 is the area of S 3 . An interesting feature here is that the finite mass black hole solutions have 
two free parameters in the event horizon initial data, which were taken to be V(rh) and V(rh)- As a result, 
and in contrast to the case of globally regular solutions, the two electric charges q and q are independent for 
black holes. This leads to a much richer parameter space of solutions. 

Our numerical results provide evidence for the existence of finite mass black hole solutions of the EYMCS 
system with a set of four notrivial gauge functions (i.e. for the full group SO(6)). The profile of a generic 
black hole corresponding to k = 0.2 and r^ — 0.5 is presented in Figure 6 (left) for q — 0.5, q = —1.5. 

However, given the large number of free input parameters, we did not attempt a systematic study of these 
solutions, concentrating instead on the simpler case of the SU(2) x U(l) truncated model. An interesting 
feature here is that solutions with AdS asymptotics exist only for the limited interval n m in < k < K max of 
k. The limits of this interval depend on the values of Th and K. Clearly, the range of the electric charge of 
these solutions is also bounded. These features are illustrated in Figure 7 (left), for black hole solutions with 
a fixed horizon radius r^ = 1. One can see that as k — > n m in the value at the horizon of metric function 
er(r) tends to zero, as does also the Hawking temperature, while other quantities stay finite. The behaviour 
of solutions for large k is less clear, the accuracy decreasing with n. However, the numerical results seem to 
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Figure 6: Left: The profiles of a typical non extremal black hole solution of the SO(6) model is presented as a 
function of the radial coordinate r. Right: An extremal black hole solution for our particular truncation of the full 
model. 



indicate that this branch ends in a critical solution with cr(r^) close to one and a finite nonzero value of Th- 
Unfortunately, the study of solutions for k — > n m ax is a difficult task and the general picture may be much 
more complicated. For example, we have noticed the existence there of a secondary branch of solutions, 
which are close to a finite mass extremal configuration with nontrivial gauge fields. A systematic study of 
these aspects would require a different parametrisation of the metric line element than ([6]), and is beyond 
the scope of this work. 

One can also keep k and K fixed and vary the value of the event horizon radius. As noted above, 
choosing the values of k,K fixes also the electric charge, i.e., these black holes are in a canonical ensemble. 
Our numerics indicate that for any K, the value of the gauge field potential at the horizon decreases with 
r/j. For large enough values of the event horizon radius, the solutions become essentially RNAdS black holes, 
with w(r) being close to the value —1 everywhere, with the non-Abelian magnetic field vanishing. 

However, the picture for small enough values of the horizon radius depends crucially on the value of 
the integration constant K in the U-equation (|29|) (i.e. on the Abelian electric charge). Starting with the 
special value K — 8k, we plot in Figure 8 a number of relevant features of the solutions, for three different 
values of the CS coefficient k. One can see that these non-Abelian black holes behave in a similar way to 
the vacuum AdS solutions. The small black holes are thermally unstable, the entropy being a decreasing 
function of the Hawking temperature. They become thermally stable for large enough values of the event 
horizon radius. In the limit r^ — > 0, these black holes approach the set of globally regular particle-like 
solutions with V(r) = w(r) = 0, discussed above . 

The picture is very different when choosing instead K 7^ 8k (see Figure 7 (right)). As an interesting 
new feature, here an extremal black hole solution is approached for a critical value of r^. In this case the 
horizon is degenerate (i.e., N(r) has a double root: N(rh) — TV'(r^) = 0) and the near horizon geometry is 
AdS 2 x S 3 . 

As r — > rh one finds the approximate form of the solution in the near horizon region: 

N(r)=N 2 (r-r h f + 0(r-r h f, <r(r) = a h - ^l( r _ rh) + (r - r h ) 2 , (32) 
w(r)=w h +w 1 (r-r h ) + 0(r-r h ) 2 , V(r) = V h - ^^(r - r h ) + 0(r - r h f , (33) 
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Figure 7: The mass parameter Mo, the value at the horizon of the metric function cr(r) and the magnetic gauge 
potential w(r), the electrostatic potential $ and the Hawking temperature are plot as a functions of the CS coupling 
constant k (left) and of the event horizon radius ru (right). 



where 



N 2 = 



Z2n 2 r h £ 2 w h {l-w 2 ) 
-3a 2 r 2 £ 2 w 2 + 384kH 2 (1 - w 2 ) 2 - 4« 2 (24r£ - 4r 2 £ 2 (w 2 - 2) + 3a 2 £ 2 (l - w 2 ) 2 ) 
244 + 8r 2 £ 2 - 3a 2 £ 2 (l - w 2 ) 2 ) 



1 ~ o^,2^2»2,„2 , OQA-Atflf-l _ „„2\2 _ A^fOA^ _ A„2p2( nn 2 _ n\ _i_ ^2i2/l _ „„2\2N ^ / 



2r£P 



The parameters r% and Wh in the above relation are solutions of the equations 

rfl 12« 2 (1 - wl) 2 + r 2 (16 - *^£) = 0, (35) 



32rt 



£ 2 

2Kk + w h (r 2 h + 8k 2 (wI - 3)) = 0. 

Recalling that K — 2(q — 4k), it follows that all event horizon boundary data (except cr{rh)) are fixed by the 
K, £ and the electric charge q. (Note the analogy with the extremal Abelian solution case.) This extremal 
solution differs from the RNAdS one, presenting no n- Abelian magnetic hair and a nontrivial metric function 
a(r) (see Figure 6 (right)). 

As expected, the near horizon structure of the extremal solutions can be extended to a full AdS2 x S 3 
solution of the field equations. This configuration has a line element 



dr o 9 / Air^ , 9 
i_ „2jr>2 (-[ 1 \j±2 

Y _ Air 2 



ds = aT^ + r o^3 - (1 - -7T-W, (36) 

6 

and the matter fields 



w = w a , V(r) = V -^-(r-r h ). (37) 

For w 2 , ^ 1, this is a non- Abelian solution, with the gauge field living on the three-sphere. The parameters 
A, wq and the radius tq of the S 3 are constrained by the relation 

3 nG r 2 w 2 + l2K 2 (w 2 -l) 2 

A - — 5 2~4 ' y 6 °) 

the value of the AdS2 cosmological constant Ai in (j3"6")l being 

A 1=6 (A-^- ^ +4 fti-«8) a y (39) 

\ ri e z K z r^ J 
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Figure 8: The temperature-entropy diagram (left) and the value at the horizon of the metric function a(r) and of 
the electrostatic potential $ (right) are ploted as functions of the event horizon radius for three different values 
of the Chern-Simons coupling constant k. These configurations have an integration constant K = 8k, and approach 
the globally regular particle-like solutions as m — > 0. 



4 Global charges 

4.1 The mass and boundary stress tensor 

The action and mass of these AdSs non-Abelian configurations is computed by using a boundary countert- 
erm prescription. As found in |20) . the following counterterms are sufficient to cancel divergences in five 
dimensions, for Schwarzschild-AdS black hole solution (in this Section we restore the 8ttG and e factors in 
the expressions): 





3 




f d i xv^h 






IBM 


T 





(40) 



with R the Ricci scalar for the boundary metric h. 

(Note also that, as usual, to ensure well-defined Euler-Lagrange field equations, one adds to the action (JTJ, 
the Gibbons- Hawking surface term [23] I SUI f — — Jqm d 4 x^—hK, where K is the trace of the extrinsic 
curvature for the boundary DM.) However, in the presence of matter fields, additional counterterms may 
be needed to regulate the action [35], which is also the case for the generic non-Abelian solutions discussed 
in the previous Section^]. 

This divergence is cancelled by a supplementary counterterm of the form (with a, b boundary indices): 

ll M = - \og(-) [ d+xV^h^ Tr{F ab F" b } . (41) 
« JdM 2e 

Note that this term is identically zero for the solutions with w 2 + w 2 — >• 1, Vw — Vw —> 0. 

Using these counterterms and the Gibbons-Hawking boundary term, one can construct a divergence-free 
boundary stress tensor T a b 

1 3 £ 2£ r 1 

T ah = 7TT;(K ab - Kh ab - -h ab + -E ab ) - - log(-) Tr{F ac F bd h cd - -h ab F cd F cd } , (42) 
»7rG I 2 e z £ 4 

11 The geometric counterterm H40H regularises also the action of the RNAdSs black hole solution. However, this does not 
hold for any d = 5 solutions of the Einstein-Maxwell-A system. An interesting example here are the AdS black strings with a 
magnetic (7(1) field 1361 , in which case one has to consider an additional matter counterterm on the form 141t . 
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where E a f, and K are the Einstein tensor and the trace of the extrinsic curvature K a t for the induced 
metric of the boundary, respectively. In this approach, the mass M of the solutions is the conserved charge 
associated with the Killing vector d/dt 20 . A straightforward computation leads to the following simple 
result for the mass of the generic EYMCS solutions: 



M 



W3M0 
16ttG 



M r , with M r 



W 3 £ 2 
64ttG 



(43) 



For the case of black hole solutions of the truncated SU(2) x U(l) model, we have found that M coincides 
within the numerical accuracy with the mass computed from the first law of thermodynamics, up to the 
constant term M c which is usually interpreted as the mass of the pure global AdSs. 

From the AdS/CFT correspondence, we expect the non-Abelian hairy black holes to be described by 
some thermal states in a dual theory formulated in a metric background given by 



jabdx a dx b 



-dt 1 



£ 2 (dip 2 + sin 2 iP(d0 2 + sin 2 8dp> 2 )), 



(44) 



where ip, 8, ip are the usual polar angles parametrizing S 3 . 

The matter fields in the dual CFT would interact with a background non-Abelian field, whose expression, 
as read from (J5J, (H3) is 



A t = 

A^ = 

Ag = 

A, n = 



- (Vo [sin ip sin 8 (£i 6 cos ip + £26 sin ip) + sin ip cos 8 £ 36 + cos ^£46] — Vo £56 
e 

1 

e 



— ((1 + wo) [sin#(Xi4 cos ip + £24 sin 93) + cos 8 £34] 
e 

+wo [cos V , sin6'(Ei5 cost/? + £25 simp) + cos ip cos #£35 — sin ^£45] ), (45) 

— ((1 + wq) sin-0 [ sinf/;(£i3 cos <p + £23 sin tp) — sin ip cos 8 £14 cos <p 
e 

+ cos ip cos 8 £24 sin tp + cos ip sin 8 £34] + wq sin tp [cos #(£15 cos (p + £25 sin ip) — sin 8 £35] ) , 

— ( — (1 + wo) sin ip \ sin ip sin 2 8 £12 + sin ip sin 8 cos #(£13 sin ip — £23 cos p) 
+ cos tp sin 6*(£i4 sin <p — £24 cos <p)] — wq sin ip sin #(£15 sin <p — £25 cos tp)) , 



We note that this is still fully an SO(6) gauge field. 

The expectation value < t£ > of the dual CFT stress tensor can be calculated using the relation [37] 



-77 Qb < r bc >= lim 



-hh ab T. 



be- 



Employing also (|42|) . we find the finite and covariantly conserved stress tensor (with x 1 = ip, x 2 = u. r 



(46) 
3 _ 



p, x^ 



8itG < r b a >= — (^r + - 



U K I 2 



/ 1 \ 

10 

10 

\ -3 J 



4irG((w 2 +w 2 - l) 2 + e 2 (V w - VqWq) 2 ) 
e 2 P 



/ 1 \ 

10 

10 

\ j 



Different e.g. from the case of Reissner-Nordstrom-AdS Abelian solutions, this stress tensor has a nonvanish- 
ing trace. Moreover, for the physically relevant case of solutions with a finite electric charge (i.e. \cp X \\ ~~ ^ 



asymptotically) one finds that < >= Aym = — 3(wq 



1) /(2re ). This agrees with the general 



results [35], [35], [3S] on the trace anomaly in the presence of an external gauge field, Aym — 
coefficient 1Z being related to the charges of the fundamental constituent fields in the dual CFT 



TlF^y the 



(47) 



4.2 Electric charge(s) 

For the solutions with \<p x \\ ~^ (*- e - (Vo^o — Vow )|r^oo = 0) the coefficients of the 1/r 2 terms in 
the asymptotic expansion of the electric potentials are finite. Thus, from the Gauss flux theorem one can 



17 



formally define the electric charge 

4tt 2 



e = <t> dS k V^F kt = — (qS 5 6 + gSia) , (48) 



DC 



(? 



which is clealy not gauge invariant. (This is a generic problem for the definition of the non-Abelian charges 
in the absence of a Higgs field, see e.g. |40) ) . 

Perhaps a more proper definition can be given following the reasoning in [41] . In this approach one starts 
by evaluating the quantity Tr-fi^F 4 *}, 



^Tr{F lt i^} - y/=$I!r {DiA t F } = d % (^Tr {A t F lt }) - ^Tr{A t A'F } ■ 
Using the Gauss' law equation, we find that the contribution of the electric field to the total mass then is 

E e = ~^J dSkV^gx-DrX + W^ J dr(\$\ 2 -l)x-D r $g tt (49) 

Subject to the truncations this expression simplifies; the integral in the second term must be evaluated using 
the numerical solution, while the surface integral in the first term can be evaluate using only the asymptotic 
values of the functions. 

In the absence of the Chern-Simons term, i.e., when k = 0, the contribution of the electric field to the 
total mass can be written as 

E e = -<f dS k ^Tr{A t F kt } = Q E * , 

J OO 

where 



* = VtHAM = ^ + v*, Q E = — nq + m , (50) 

v e ^J V 2 + V 2 

are the electrostatic potential and the electric charge, respectively. However, one can extend this definition 
of $ and Qe to solutions of the EYMCS system. This applies to both generic and special configurations 
(note that $ = $, Qe = 47r 2 (sina q + cos a q)/e for finite mass solutions). 

For finite mass solutions, following |21) . one can also define an effective non-Abelian charge Q e ff by the 
asymptotic behaviour of the metric function N(r), which, to order 1/r 4 is similar to that of the RNAdS 
solution: 



^>=M---^ + %+.... (.VI, 



a^{g 2 + q 2 )P + i{wl+w 2 2 )) 

Qeff = J ■ (52) 

Concerning a definition of a "magnetic" flux, the only natural quantity we have at our disposal for this 
purpose is the Chern-Pontryagin density, which we know is the leading and sole contributing term to the 
topological charge ("magnetic" flux) of the monopole in 4 + 1 dimensions [42] . There however the gauge 
group is 50(4) and the model features a iso-four-vector Higgs field. 

This quantity can be calculated easily for the SO(6) Ansatz employed in this paper 

e ijk l Tr {F 13 F kl } = ~ {\$\ 2 - 1) Tr (D P <f> 3 S 56 + (D r cpe) M E M4 ) , (53) 

which vanishes, i.e., the candidate for a magnetic charge for the solutions found in the present work equals 
zero identically. 
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5 Further remarks 



On general grounds, one expects that extending the known classes of solutions of the d — 5 supergravity to 
a non-Abelian gauge group would lead to a variety of new physical effects. 

This work has been aimed as a first step towards constructing the non-Abelian solutions of the maximal 
d = 5 gauged supergravity. Restricting to the simplest case of static, spherically symmetric solutions, we 
have proposed a suitable Ansatz for the gauge fields and presented numerical evidence for the existence of 
both particle-like and black hole solutions. Our systematic description of the black holes is restricted to a 
certain truncation of the full 5*0(6) model, with the sole purpose of rendering the numerics practicable. In 
this limited context, we have also found extremal black holes. 

As a consequence of the presence of a negative cosmological constant in the model, we have recovered 
the qualitative properties of the solutions to the usual EYM- A model in 3 + 1 dimensions pQ , [2] . Notably, 
some of our solutions to which we have referred as generic, are characterised by arbitrary asymptotic values 
of the potentials parametrising the gauge field. Also, these solutions share another property with those of 
the 3 + 1 dimensional EYM model, namely that the shooting parameters involved take on a continuum of 
values. Unlike the latter however, their masses turn out to be divergent in our 4 + 1 dimensional case. This 
is expected on the basis of the Derrick-type scaling argument. However, we have proposed a regularisation 
procedure for the mass of these solutions, in the context of the AdS/CFT correspondence. As far as these 
generic solutions are concerned, the presence of the Chern-Simons term makes no qualitative difference. 

Perhaps the most interesting feature of the EYMCS model is the existence of finite mass solutions. We 
have referred to these as special solutions and they contrast with the generic ones in that the shooting 
parameters involved take on a discrete set of values. The special solutions exist only when the non-Abelian 
Chern-Simons term is present. 

Concerning the physical context of our results, it is in order to make several remarks on the issue of the 
more general solutions of the d — 5, M = 8 gauged supergravity. This model contains in addition twenty 
scalars, which are represented by a symmetric unimodular tensor. These scalars have a nontrivial potential 
approaching a constant negative value at infinity which fixes the value of the effective cosmological constant. 
No obvious consistent truncation of this sector seems to exist for a gauge group SO(6) and strictly speaking- 
one should work with the full set of scalars. In principle, at least when this general model is subjected to 
spherical symmetry, the one dimensional subsystem resulting from the aplication of the Ansatz here can be 
studied using the same methods as in this paper, i.e., the solutions can be constructed by solving a boundary 
value problem. One can in that case find the approximate expressions at the origin or event horizon and at 
infinity. The only obstacle to this task that we see at this moment is the huge complexity of the ensuing 
equations. Based on the results in this paper, we expect the parameter space of the full SO(6) solutions of 
the d — 5, J\f — 8 model to be very rich. Inclusion of the scalar sector will lead to many new free parameters 
in the asymptotics and will make any attempt to classify the solutions difficult (involving numerous different 
ways of approaching a constant negative value at infinity for the scalar potential). 

The generic non-Abelian solutions on the other hand will always present a nonvanishing magnetic gauge 
field on the boundary which appears as a background for the dual theory. (This feature is independent of 
the presence or absence of scalars.) Thus the expectation value of the dual CFT stress tensor will contain 
a part which is similar to (|47|) . Configurations with vanishing non-Abelian magnetic field on the boundary, 
and finite mass, should exist as well, being supported by the CS term. Also, similar to the case of four 
dimensional EYM non-Abelian solutions in [6] , the existence of both spherically symmetric globally regular 
and hairy black hole solutions with the same set of data at infinity raises the question as to how the dual 
CFT is able to distinguish between these different bulk configurations. 

Concerning future developments of this work, we expect a much richer structure of the non-Abelian so- 
lutions to be found when relaxing the spacetime symmetries. For example, one may envisage the existence 
of asymptotically AdSs solutions which are static and non-spherically symmetric, generalising the configu- 
rations in [43] ■ Particularly interesting would be to approach the issue of EYMCS rotating solutions. To 
our knowledge, the only d > 4 rotating solutions with non-Abelian fields known so far in the literature are 
the d — 5 EYM- SU (2) black holes in [44]. However, the mass of these solutions diverges logarithmically. It 
is likely that the inclusion of a CS term will lead to finite mass solutions also in that case. 
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Another very natural direction to be explored is the case of zero cosmological constant, which would 
be outside the context of M = 8 supergravity but would nonetheless be technically very interesting. For 
example, this would afford a comparison with the corresponding 3 + 1 dimensional asymptotically flat EYM 
solutions in |21) . which involve only a discrete spectrum of the shooting parameters. 
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